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ABSTRACT 

We present spectral line profiles for the 557 GHz 11,0 ~ ► lo.i ground-state rotational 
transition of ortho-H 2 16 0 for 18 galactic star formation regions observed by SWAS. 

Water is unambiguously detected in every source. The line profiles exhibit a wide 
variety of shapes, including single-peaked spectra and self-reversed profiles. We 
interpret these profiles using a Monte Carlo code to model the radiative transport. The 
observed variations in the line profiles can be explained by variations in the relative 
strengths of the bulk flow and small-scale turbulent motions within the clouds. Bulk 
flow (infall, outflow) must be present in some cloud cores, and in certain cases this 
bulk flow dominates the turbulent motions. 

Subject headings: ISMrmolecules - radio lines:ISM - stars:formation - ISMxlouds 


1. Introduction 


Water in interstellar space is of great interest for a variety of reasons, among them that it 
is widely believed to be a significant participant in the chemistry of molecular clouds (Bergin, 
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Neufeld, Sz Melnick 1998) and a major coolant during the later stages of star formation (Ceccarelli, 
Hollenbach, Sz Tielens 1996). Unfortunately, most H 2 0 transitions are obscured by the Earth’s 
atmosphere. Recently the Infrared Space Observatory (ISO) has detected high-energy H 2 0 
transitions towards numerous dense cores (Liseau et al. 1996, Helmich, et al. 1996, van Dishoeck 
& Helmich, 1996, and Cernichaxo et al. 1997), but these transitions tend to be predominantly 
sensitive to hot gas surrounding young stars (van Dishoeck et al. 1999). 

Observations by the recently launched Submillimeter Wave Astronomy Satellite (SWAS) of 
the 556.936 GHz 1 1(0 -► l 0 ,i ground state transition of H 2 0 have opened another window on 
water in molecular clouds. This transition, with an upper state energy only 27 K above the 
ground state, is predicted to be a good probe of the excitation conditions within the extended 
colder gas surrounding young stars (cf. Liseau & Oloffson 1999). In addition, the high optical 
depths typically encountered in this line give rise to strong radiative trapping. This combination 
of traits can lead to significant self-absorption and complex H 2 0 profiles which contain important 
information about cloud structure, i.e., the density, temperature, and velocity structure, and even 
cloud geometry. However, very optically thick lines often fail to exhibit self reversals even in the 
presence of excitation gradients. This is a problem of long standing in molecular astronomy, which 
has been addressed by models that vary the relative contributions of local turbulent motions and 
the bulk velocity field (Goldreich & Kwan 1974, Scoville & Solomon 1974, Martin, Sanders, Sz Hills 
1984, Kwan & Sanders 1986, Juvela 1997). These models reduce the magnitude of self-absorption 
in the optically thick lines. 

In this Letter we present SWAS observations of the 557 GHz 11,0 lo,i transition of 
ortho-H 2 0 toward 18 cloud cores. We use spherically symmetric Monte Carlo radiative transfer 
models to examine the cloud structures that produce these complex spectra. 


2 . Observations 

SWAS has observed numerous dense cores in molecular clouds. The sources presented in 
this work tend to be massive cloud cores with active star formation and, often, embedded stellar 
clusters. For a discussion of H 2 0 emission in less massive cores, see Snell et al. (2000a) and Ashby 
et al. (2000). SWAS observations of Orion BN/KL are presented in Melnick et al. (2000b), while 
extended H 2 0 emission in OMC-1 is discussed in Snell et al. (2000b). 

The data were collected in the standard nodding mode (Melnick et al. 2000a) in segments 
typically 30 minutes long during 1999 January - October. SWAS has a main beam efficiency of 
0.9, an average beam width of 3. 7 8 (FWHM), and a velocity resolution of 0.6 km s 1 (Melnick et al. 
2000a). Source positions and distances are given in Table 1. Figure 1 presents the SWAS spectra. 
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3. Results 

SWAS has detected H 2 0 toward every object in Table 1. H 2 0 emission appears to be 
ubiquitous in massive cloud cores. The majority of these sources exhibit double-peaked spectra 
with deep troughs. Most of these troughs coincide closely with the systemic velocity inferred 
from the simultaneously observed ^CO J = 5 — > 4 line, indicating strong self- absorption in the 
H 2 0 line. The others, (e.g., DR 21 and NGC 63341), are due to absorption by foreground material. 
Several sources which might superficially appear to be dominated by a single component (e.g., 

G 265.1+1.5, NGC 2071, and NGC 7538) differ from the ‘truly’ single-peaked sources in that their 
H 2 0 emission peaks are significantly offset from the systemic velocity; we infer they are likewise 
self-reversed. The spectra presented in this Letter fall roughly evenly into three general categories. 
Five spectra are single-peaked; six are self-reversed with broad lines (FWHM> 5 km s ), and 
seven are self-reversed with relatively narrow lines. 

This behavior appears to be reflected by 4 * * * * * * * 12 CO line profiles in the literature. Sources with 
self-reversed H 2 0 lines tend to exhibit self-reversed 12 CO (2 1) profiles. G 265.1+1.5, NGC 2071, 

and NGC 7538, which superficially appear to be single-peaked, likewise exhibit such reversals 
(Brand et al. 1984, Phillips et al. 1981, and Kameya et al. 1989). Conversely, our single-peaked 
sources (G 291.3-0.7, G 322.2+0.6, S140) are single-peaked in 12 CO (2 1) (Brand et al. 1984, 

Phillips et al. 1988). M17SW lacks a published 12 CO (2 -> 1) spectrum, but its 3 ^ 2 and 
1 — v 0 spectra are both single-peaked (Martin, Sanders, k Hills 1984). Thus it appears that the 
physical conditions affecting the * 2 CO line influence the H 2 0 line in an analogous manner. This is 
somewhat surprising since 12 CO and H 2 0 have such different critical densities. 

Significant broad (10-15 km s -1 ) H 2 0 emission is detected in most sources, centered on the 
systemic cloud velocities. These broad wing features are presumably from H 2 0 entrained in 
known high-velocity outflows present in the 3^3 x 4f 5 SWAS beam, e.g., p Oph A, Ceph A, DR 21, 
NGC 1333, NGC 2071, NGC 63341, NGC 7538, and S140. See Neufeld et al. (2000) for a discussion 
of SWAS spectra of molecular outflows. 


4. Monte Carlo Models of H 2 0 Emission 

The variety of H 2 0 line profiles shown in Figure 1 suggests that there are differing physical 

conditions within their clouds of origin. Since many of the clouds have quite similar density (Snell 

et al. 1984, Mundy et al. 1986, Bergin, Snell, k Goldsmith 1996, and Snell et al. 2000a), and 

temperature (Bergin et al. 1994, Churchwell k Hollis 1983) structures, velocity structure likely 

plays an important role in line formation (cf. Kwan k Sanders 1986). 

To test this idea we have developed a Monte Carlo radiative transfer code to examine 

how physical structure affects line formation for the 557 GHz transition. The Monte Carlo 

treatment provides a careful global accounting of photon absorption and reemission that is missing 
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from simpler (but computationally faster) models such as the Large Velocity Gradient (LVG) 
approximation. Our code is a modified version of the original Bernes (1979) code, to which we 
added a treatment of continuum photons from dust mixed with the gas. Our code has been 
extensively tested by comparing its results to those obtained using escape probability methods 
in the LVG limit and to those generated by the Accelerated Lambda Iteration (ALI) code of 
J. Zmuidzinas in the static limit. The code uses a physical model consisting of a spherically 
symmetric cloud subdivided into an arbitrary number of shells. Physical conditions may vary 
within a cloud but are held constant within each shell. The ability to vary physical parameters 
arbitrarily affords great flexibility. 

This code allows us to determine how different contributions from turbulent and outflow 
motions affect the simulated line profile. Our approach was to construct a family of simple 
constant temperature and density models (T = 30 K, n Hz = 10 5 cm" 3 ) of spherical outflow and 
infall. For the outflow models (Figure 2), the velocity, v, is zero at the cloud center and increases 
linearly to vrmax at the cloud edge. To create infall models (Figure 3) we subtracted vrmax 
from the outflow model velocity fields, which yields a velocity equal to zero at the cloud edge 
and decreasing toward the cloud center. We varied the relative strengths of the local turbulent 
velocity width vturb and the bulk velocity field while holding the sum of the two contributions 
constant, i.e., V 0 = (v% URB + (2 vrmax ) 2 ) 112 for outflow and V 0 = {v$ URB + (2 v C en ) 2 ) 1/2 for ^ 
infall, where vcen is the velocity at the cloud center. We chose a total width of Vo = 4.0 km s , 
and set RMAX = 0.2 pc, these values being typical for the sources in Figure 1. We assume a 
power law of the form k = ^(A/Ao)" 1 ' 5 cmV 1 for the opacity per gram of dust and gas, where 
A 0 = 1300 /im and k 0 = 5x 10 -3 cm 2 g _1 (Goldsmith et al. 1997, Hildebrand 1983). Within each 
shell we interpolate Phillips et al. (1996) collision rates for ortho- and para-H 2 to derive rates 
appropriate for the given temperature of 30 K. We assume a ratio of ortho- to para-H 2 of 0.1, and 
a volume filling factor of unity. We modeled five energy levels in the ortho-H 2 0 molecule, i.e., all 
levels up to T~ 200 K. 

Though these models are simplistic compared to more physical treatments (e.g., Shu (1977) 
infall), they are nonetheless effective in illustrating the effects of bulk flow on the H 2 0 spectra. 
The model spectra exhibit single-peaked lines when bulk flows dominate. Self-absorption results 
when turbulent motions dominate. When both components of the velocity field are comparable, 
an obvious asymmetry is apparent. For the model spectra in Figure 2, the blue peak is weaker, as 
expected, for outflow models. The opposite asymmetry occurs for the infall models. Though bulk 
flows appear to dominate in most cases, the model spectra presented here exhibit the full range of 
profiles seen in the SWAS data. For example, NGC2024 has a spectrum that may be compared to 
the vturb/Zvrmax = 0.75 model, i.e., a model in which both contributions to the velocity field 
are comparable. On the other hand, any of the models in which bulk flow is dominant appears 
to agree better with the 557 GHz H 2 0 spectrum of M17SW. Finally, infall models with a small 
contribution from turbulent motions resemble the spectra of G267.9-1.1 and NGC 63341. Detailed 
models of S 140 and p Oph A are discussed in Ashby et al. (2000). 
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5. Summary 

The SWAS spectra presented here demonstrate that 557 GHz H 2 O emission is common in 
massive cloud cores. The line profiles vary greatly from source to source, but a common feature 
is a broad component attributable to known outflows within the S"WAS beam. In addition, the 
spectra are typically either single-peaked, self-absorbed with broad components, or self-absorbed 
with relatively narrow components. CO line profiles may be useful as predictors of the ortho-H 2 0 
line profiles. The profiles’ variety can be attributed to differences in the relative contributions of 
the bulk flow and the turbulent motions to the overall velocity field; bulk flows may in some cases 
dominate the turbulent motions for the 557 GHz water line. 

The authors gratefully acknowledge helpful discussions with P. Myers and M. Crosas. J. 
Zmuidzinas generously made his ALI code available. This work was supported by NASA Grants 
NAS5-30702, NAG5-3542, and NASA’s Long Term Space Astrophysics Program. 
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Table 1. 


SOURCE 
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Targets Observed by SWAS 


POSITION (J2000) DISTANCE 


6 (pc) 


2:27:02.8 61:52:21 2200 

3:29:10.5 31:13:35 350 

5:35:14.5 -5:25:49 415 

5:41:44.5 -1:55:35 415 

5:47:04.1 0:21:43 390 

6:41:03.9 9:34:39 800 

8:59:12.0 -47:29:04 1700 

8:59:23.3 -43:46:14 1000 

9:01:54.3 -47:43:59 800 

11:11:37.7 -61:19:49 2400 

15:18:35.2 -56:38:25 3800 

16:26:23.4 -24:23:02 160 

17:20:53.0 -35:44:57 1740 

18:20:22.1 -16:12:37 2200 

20:39:00.9 42:19:38 3000 

22:19:17.1 63:18:46 910 

22:56:17.9 62:01:50 730 

23:13:47.6 61:26:54 
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i*S 




pig i — SWAS spectra of the 557 GHz transition of ortho-H20. Dotted vertical lines mark the 
fitted line center velocity of the simultaneously observed 550.926 GHz 13 CO J = 5 -> 4 transition. 
OMC-l(0,-3.2) was observed at a =5:35:14.5, <5 =-5:25:49 (J2000), 3f2 south of the BN/KL object. 
A first-order baseline has been subtracted from the data. Double-peaked spectra have been placed 
in the right-hand (relatively broad lines) and center (relatively narrow lines) columns, and the 
left-hand column contains single-peaked spectra (but see the discussion of G265.1+1.5 in the text). 
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VELOCITY/V 0 


Fig. 2. — Monte Carlo model 557 GHz H 2 O spectra for spherical outflow in clouds with constant 
temperature and density, in which the relative contributions of the turbulent velocity and the linear 
spherical outflow velocity field to the total velocity V 0 are varied as described in the text. 






